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Abstract 



We study the double charm decays of B c meson, by employing the perturbative QCD approach 

based on kx factorization. In this approach, we include the non-factorizable emission diagrams 
<D ■ 
,-C and W annihilation diagrams, which are neglected in the previous naive factorization approach. 

pg ■ The former are important in the color-suppressed modes; while the latter are important in most 

> ■ 

£T) • B c decay channels due to the large Cabibbo-Kobayashi-Maskawa matrix elements. We make 

o. 

comparison with those previous naive factorization results for the branching ratios and also 



give out the theoretical errors that previously missed. We predict the transverse polarization 



fractions of B c — > D*,\ D*° , D?~t D*° decays for the first time. A large transverse polarization 
contribution that can reach 50% ~ 60% is predicted in some of the B c meson decays. 

j_l i PACS numbers: 13.25.Hw, 12.38.Bx, 14.40.Nd 



* Electronic address: lucd@ihep.acxn] 



I. INTRODUCTION 

Since the B c meson is the lowest bound state of two different heavy quarks with open 
flavor, it is stable against strong and electromagnetic annihilation processes. The B c 
meson therefore decays weakly. Furthermore, the B c meson has a sufficiently large mass, 
thus each of the two heavy quarks can decay individually. It has rich decay channels, and 
provides a very good place to study nonleptonic weak decays of heavy mesons, to test the 



Q. 



The current running LHC 



standard model and to search for any new physics signals 

collider will produce much more B c mesons than ever before to make this study a bright 

future. 

Within the standard model (SM), for the double charm decays of B u ^ s mesons, there 
are penguin operator contributions as well as tree operator contributions. Thus the di- 
rect CP asymmetry may be present. However, the double charm decays of B c meson 
are pure tree decay modes, which are particularly well suited to extract the Cabibbo- 



Kobayashi-Maskawa (CKM) angles due to the absented interference from 



aenguin contri- 
6[, the decays 



butions. As was pointed out in ref. [2] and further elaborated in ref. [SHS] 
B c —¥ D^D°, DfD° are the gold-plated modes for the extraction of CKM angle 7 though 
amplitude relations because their decay widths are expected to be at the same order of 
magnitude. But this needs to be examined by faithful calculations. 

Although many investigations on the decays of B c to double-charm states have been 



carried out 



Q,Q,QJi2 



12J in the literature, there are uncontrolled large theoretical errors with 
quite different numerical results. In fact, all of these old calculations are based on naive 
factorization hypothesis, with various form factor inputs. Most of them even did not give 
any theoretical error estimates because of the non-reliability of these models. Recently, 
the theory of non-leptonic B decays has been improved quite significantly. Factorization 
has been proved in many of these decays, thus allow us to give reliable calculations of 
the hadronic B decays. It is also shown that the non-factorizable contributions and 
annihilation type contributions, which are neglected in the naive factorization approach, 



are very important in these decays 13 ] 



The perturbative QCD approach (pQCD) [l4| is one of the recently developed the- 
oretical tools based on QCD to deal with the non-leptonic B decays. Utilizing the k^ 
factorization instead of collinear factorization, this approach is free of end-point singular- 



ity. Thus the Feynman diagrams including factorizable, non-factorizable and annihilation 
type, are all calculable. Phenomenologically, the pQCD approach successfully predict the 
charmless two-body B decays [15|, LL6( . For the decays with a single heavy D meson in the 
final states (the momentum of the D meson is |mg(l — r 2 ), with r = mp/m^), it is also 
proved factorization in the soft-collinear effective theory [17] . Phenomenologically the 
pQCD approach is also demonstrated to be applicable in the leading order of the m D /m B 



expansion 18|, |l9] for this kind of decays. For the double charm decays of B c meson, the 
momentum of the final state D meson is |me c (l — 2r 2 ), which is only slightly smaller 
than that of the decays with a single D meson final state. The prove of factorization here 
is thus trivial. The pQCD approach is applicable to this kind of decays. In fact, the 
double charm decays of B u ^, s meson have been studied in the pQCD approach success- 



fully 20|, [2l| , with best agreement with experiments. In this paper, we will extend our 
study to these B c decays in the pQCD approach, in order to give predictions on branching 
ratios and polarization fractions for the experiments to test. Since this study is based on 
QCD and perturbative expansion, the theoretical error will be controllable than any of 
the model calculations. 

Our paper is organized as follows: We review the pQCD factorization approach and 
then perform the perturbative calculations for these considered decay channels in SecJTTl 
The numerical results and discussions on the observables are given in Sec lIHI The final 
section is devoted to our conclusions. Some details of related functions and the decay 
amplitudes are given in the Appendix. 

II. FRAMEWORK 

For the double charm decays of B c , only the tree operators of the standard effective 
weak Hamiltonian contribute. We can divide them into two groups: CKM favored decays 
with both emission and annihilation contributions and pure emission type decays, which 
are CKM suppressed. For the former modes, the Hamiltonian is given by: 

H eff = ^v^v^cmo^) + cmom), 

2 = b a ^(l - 7 5 )c a ® %7 M (1 - 75)9/3, (1) 



while the effective Hamiltonian of the latter modes reads 

H e ff = ^v^v^cmo^ + cmo'M}, 

0[ = 1^(1 - 75)% ® c^7 M (l - 7 5 )g a , 

0' 2 = 6 a 7 M (l - 7s)w Q ® C/^ 1 - 75)9/3, (2) 

where V^(g = d, s) are the corresponding CKM matrix elements, a, (3 are the color 
indices. C\^ are Wilson coefficients at renormalization scale \x. O12 and 0[ 2 are the 
effective four-quark operators. 

The factorization theorem allows us to factorize the decay amplitude into the convolu- 
tion of the hard subamplitude, the Wilson coefficient and the meson wave functions, all 
of which are well-defined and gauge invariant. It is expressed as 

C(t)®H(x,t)®Q(x)®exp[-s(P,b)-2 [ —J q {a s {n))), (3) 

Jl/b ^ 

where C(t) are the corresponding Wilson coefficients of effective operators defined in 
eq. ( [11121) . Since the transverse momentum of quark is kept in the pQCD approach, the 
large double logarithm In (Pb) (with P denoting the longitudinal momentum, and b the 
conjugate variable of the transverse momentum) to spoil the perturbative expansion. A 
resummation is thus needed to give a Sudakov factor exp[— s(P, b)] 22]. The term after 
Sudakov is from renormalization group running with 7 g = —q: s /tt the quark anomalous di- 
mension in axial gauge and t the factorization scale. All non-perturbative components are 
organized in the form of hadron wave functions $(x) (with x the longitudinal momentum 
fraction of valence quark inside the meson), which can be extracted from experimental 
data or other non-perturbative methods. Since the universal non-perturbative dynam- 
ics has been factored out, one can evaluate all possible Feynman diagrams for the hard 
subamplitude H (x, t) straightforwardly, which include both traditional factorizable and 
so-called "non-factorizable" contributions. Factorizable and non-factorizable annihilation 
type diagrams are also calculable without end-point singularity. 

A. Channels with both emission and annihilation contributions 

At leading order, there are eight kinds of Feynman diagrams contributing to this type 
of CKM favored decays according to eq.([T]). Here, we take the decay B c — > D + D° as 
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FIG. 1: Feynman diagrams for B c — > D + D° decays. 



an example, whose Feynman diagrams are shown in Fig{TJ The first line are the emis- 
sion type diagrams, with the first two contributing to the usual form factor; the last 
two so-called "non-factorizable" diagrams. In fact, the first two diagrams are the only 
contributions calculated in the naive factorization approach. The second line are the an- 
nihilation type diagrams, with the first two factorizable; the last two non-factorizable. 
The decay amplitude of factorizable diagrams (a) and (b) in FigJT] is 



T e = -2^C F f B f 3 nM B f dx 2 f hhdhdb 



2 <p 2 {x 2 ) exp( 



b\u 2 



B 



X 



{[-(r 2 - 2)r b + 2r 2 x 2 - x 2 ]a s (t a )h e (a e , (3 a , h, b 2 )S t {x 2 ) exp[-S ab (t a )} 
+2r 2 a s {t b )h e (a e , f3 b ,b 2 ,b 1 )S t (x 1 ) exjp[-Sab{t b )]}, (4) 

where r b = m b /MB, Tj = rrii/MB(i = 2,3) with m 2) m 3 are the masses of the recoiling 
charmed meson and the emitting charmed meson, respectively; CV = 4/3 is a color factor; 
/3 is the decay constant of the charmed meson, which emitted from the weak vertex. The 
factorization scales t a ^ b are chosen as the maximal virtuality of internal particles in the 
hard amplitude, in order to suppress the higher order corrections 23| . The function h e 
and the Sudakov factor exp[— S] are displayed in the Appendix El D meson distribution 



amplitude <fi{x) are given in Appendix O The factor S t (x) is the jet 



from the threshold resummation, whose definitions can be found in 24]. 



unction resulting 



The formula for non-factorizable emission diagrams Fig. [T] (c) and (d) contain the 



kinematics variables of all the three mesons. Its expression is: 

M e = --C F f B irM B / dx 2 dx 3 / b 2 b 3 db 2 db 3 (f) 2 (x 2 )(f) 3 (x 3 ) exp( — 2 — -) x 
{[1 - xi - x 3 - r 2 (l - x 2 )]a s (£ c )/i e (/3 c , «e, &3, 62) expf-S^^)] - 
[1 - xi - x 2 + x 3 - r 2 (l - x 2 )]a s (t d )h e ((3 d , a e , b 3 , b 2 ) exp[-S cd (t d )]}. (5) 

Generally, for charmless decays of B meson, the non-factorizable contributions of the 
emission diagrams are small due to the cancelation between Fig. H](c) and (d). While for 
double charm decays with the light meson replaced by a charmed meson, since the heavy 
c quark and the light quark is not symmetric, the non-factorizable emission diagrams 
ought to give remarkable contributions. This has been shown in the pQCD calculation 
of B — > Dn decays for a very large branching ratios of color-suppressed modes [25[ and 
proved by the B factory experiments. 

The decay amplitude of factorizable annihilation diagrams Fig. [1] (e) and (f ) involve 
only the two final states charmed meson wave functions, shown as 

/»1 /»oo 

7 a = -8C F f B TTMg / dx 2 dx 3 / b 2 b 3 db 2 db 3 (j) 2 (x 2 )(i) 3 (x 3 ) x 
Jo Jo 

{[1 - x 2 }a s (t e )h e (a a , f3 e , b 2 , b 3 ) exp[-S ef (t e )}S t (x 3 ) - 

[1 - x 3 ]a s (t f )h e (a a , (3 f , b 3 , b 2 ) exjp[-S e f(tf)]S t (x 2 )}. (6) 



For the non-factorizable annihilation diagrams Fig. [T](g) and (h), the decay amplitude is 

M a = -C F f B TiM B I dx 2 dx 3 I 6162^1^2^2(^2)03(^3) exp(- 
•j Jo Jo 

x{[xi +x 3 - 1 - r c ]a s (tg)h e (/3 g , a a ,b!,b 2 ) exp[-S gh (t g )] 



h 2 l ! 2 



-[n - x 2 }a s (t h )h e ((3 h , a a , h, b 2 ) exp[-S gh (t h )}}, (7) 

where r c = itl c /Mbi with m c the mass of c quark in B c meson. Finally, the total decay 
amplitude for B c — > D + D° can be given by 

A(B C -)• D + D°) = V; b V ud [a 2 7 e + C 2 M e + a x T a + C x M a ], (8) 

with the combinations of Wilson coefficients ai = C 2 + Ci/3 and a 2 = C\ + C 2 /3, charac- 
terizing the color favored contribution and the color-suppressed contribution in the naive 
factorization, respectively. The total decay amplitudes of B c — > D+D°, B c — > D + D*° and 



B c — y DfD*° can be obtained from eq.flH]) with the following replacement: 

A(B C -> D+D°) = V; b V us [a 2 F e + C 2 M e + a^ a + C x M a }\ D+ ^ Dt , 
A(B C ^D+D*°) = V: b V ud [a2Fe + C 2 M e + a 1 F a + C 1 M a }\ D o^ D *o, 
A(B C -)■ D^D* ) = V; b V us [a 2 F e + C 2 M e + a x 7 a + C x M a \ \ D+ -+ D f,D^D*°- (9) 

Comparing our eq. flH1l9|) with the formulas of previous naive factorization approach [4J,|5|,|7- 
llOl ]. it is easy to see that only the first term appearing in eq. (!8|9|) are calculated in the 
previous naive factorization approach. The second, third and fourth terms in these equa- 
tions, are the corresponding non-factorizable emission type contribution, factorizable and 
non-factorizable annihilation type contributions, respectively, which are all new calcula- 
tions. 

In B c — y D*,\D*° decays, the two vector mesons in the final states have the same helicity 
due to angular momentum conservation, therefore only three different polarization states, 
one longitudinal and two transverse for both vector mesons, are possible. The decay 
amplitude can be decomposed as 

A = A L + A N e T 2 ■ el + M^nVe^, (10) 

where e^, ej are the transverse polarization vectors for the two vector charmed mesons, 
respectively. A L corresponds to the contributions of longitudinal polarization; A N and 
A T corresponds to the contributions of normal and transverse polarization, respectively. 
And the total amplitudes A L,N,T have the same structures as eq.f l8|9|) . The factorization 
formulae for the longitudinal, normal and transverse polarizations are listed in Appendix 

El 

For B c — y D*AD° decays, only the longitudinal polarization of D*A meson will con- 
tribute, due to the angular momentum conservation. We can obtain their decay ampli- 
tudes from the longitudinal polarization amplitudes for the B c —y D^tD* decays with 
the replacement D*° —y D°. 



B. Channels with pure emission type decays 

There are also eight kinds of Feynman diagrams contributing to B c —y D\*! £)W° 
decays according to eq.([2]), but all are emission type. Taking the decay B c — y D + D° as 
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FIG. 2: Color-suppressed emission diagrams contributing to the B c — > D + D° decays. 
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FIG. 3: Color-favored emission diagrams contributing to the B c — > D + D° decays. 



an example, Fig. |2]are the color-suppressed emission diagrams while Fig. [3] are the color- 
favored emission diagrams. We mark the subscript 2 and 3 to denote the contributions 
from Fig. [2] and Fig. EJ respectively. The decay amplitude of factorization emission 
diagrams J 7 ^, coming from Fig. [2] (a,b), is similar to eq.(jl]), but with the replacement 
D° — > D°. While the decay amplitude of non- factorization emission diagram Jvl e2 , coming 
from Fig. [2] (c,d), is different from eq.Q, since the heavy c quark and the light anti-quark 
are not symmetric. The expression of the non-factorizable emission diagram is 



h 2 , i 2 



Me2 = --C F f B irM B / dx 2 dx 3 / b 2 b 3 db 2 db 3 (() 2 (x 2 )(f> 3 (x 3 ) exp( 
<J Jo Jo 

x{[2 -xi-x 2 -x 3 -r 2 (l - x 2 )]a s (t c )h e ((3 c ,a e ,b 3 ,b 2 ) exp[-S cd (t c )] - 

[x 3 - xi - r 2 (l - x 2 )]a s (t d )h e (/3 d , a e , b 3 , b 2 ) exp[-S cd (t d )}}. (11) 

By exchanging the two final states charmed mesons in Fig. |2} one can obtain the corre- 
sponding decay amplitudes formulae T e3 and M. e3 for Fig. [31 The total decay amplitude 
of B c — > D + D° decay can be written as 



A(B C -> D + D°) = V: h V cd \a 2 T e2 + C 2 M e2 + a x T e3 + C X M 



e3 



(12) 



If the final recoiling meson is the vector D* meson, the decay amplitudes of factorization 
emission diagrams and non-factorization emission diagrams are given as 



-2^C F f B f 3 7rM B J dx 2 J b 1 b 2 db 1 db 2 



J 7 ,, 2 x l^C r fuf^Ml; I (I.ro I h i b, l lb ll ll,y, 2 (x 2 )exp(-^: 



x{[-(r 2 - 2)r b + 2r 2 x 2 - x 2 ]a s {t a )h e (a e , (3 a , b x , b 2 )S t (x 2 ) exp[-S ab (t a )] 
+r 2 2 a s {t b )h e (a e , (3 b , b 2 , bi)S t (xx) exp[-S ab (t b )]}, (13) 

8 f 1 f°° b 2 cu 2 

M* e2 = --C F f B TtM% I dx 2 dx 3 / b 2 b 3 db 2 db 3 <fi 2 (x 2 )(j) 3 (x 3 )exp( — 2 — -) 

-J Jo Jo z 

x{[2 - xi - x 2 - x 3 -r 2 (l - x 2 )]a s (t c )h e ((3 c ,a e ,b 3 ,b 2 ) exp[-S cd (t c )] - 



[x 3 - x x + r 2 (l - x 2 )]a s (t d )h e ((3 d , a e , b 3 , b 2 ) exp[-S cd (t d )}}. 



The total decay amplitudes for other pure emission type decays are then 

A(B C -)■ DtD°) = V: b V cs \a 2 T e2 + C 2 M e2 + a x T e3 + C x M e3 \ 
A(B C ^D + D*°) = V: b V cd [a 2 F e2 + C 2 M e2 + a 1 F: 3 + C 1 M* e3 } J 
A(B C ^D*+D°) = V: b V cd \a 2 r e2 + C 2 Ml 2 + a x T e3 + C x M e3 l 



(14) 



A(B C -> DfD*») = V: b V cs [a 2 F e2 + C 2 M e2 + a x r e3 + C x Mi 



e3h 



A(B C -> D* s + D°) = V: b V cs [a 2 r e2 + C 2 M* e2 + a x F e3 + CxMe^. 



(15) 



The B c — > DJtD* decays have a similar situation to B c — > D*AD*°, their factorization 
formulae are also listed in Appendix jAl 



III. NUMERICAL RESULTS 

In this section, we summarize the numerical results and analysis in the double charm 
decays of the B c meson. Some input parameters needed in the pQCD calculation are 
listed in Table [D 

A. The Form Factors 

The diagrams (a) and (b) in FigJTJor Figj3]give the contribution for B c — y DYJ. tran- 
sition form factor at q 2 = point. Our predictions of the form factors are collected in 

9 



TABLE I: Parameters we used in numerical calculation 26J 



Mass(GeV) 


M w = 80.399 


M Bc = 6.277 m b = 4.2 m c = 1.27 


CKM 


\V uh \ = (3.47lg;i|) x 10- 3 

IT/ 1 — n Q7-Ut;+0-00015 

\v cs \ - u.y^4b_ 000016 


IT/ 1 _ n Q7/I9S+0-00015 it/ 1 _ n oor: Q+0.0007 
1 Vud\ - U.9^42»_ 00015 | V us \ - U.ZZb6_ 0007 

IT/J - n 2252+ ' 0007 IT/ J - o 041 n+o- 0011 

I ''cdl — U.zzjz_ 0.0007 1 ^ cfe l — U-LKtlu -0.0007 


Decay constants (Me V) 


/fie = 


= 489 / D = 206.7 ±8.9 f Da = 257.5 ± 6.1 


Lifetime 


TBc = 0.453 x 10" 12 s 



TABLE II: The form factors for B c 



D(V at q 2 = evaluated in the pQCD approach. The 



uncertainties are from the hadronic parameters. For comparison, we also cite the theoretical 
estimates of other models. 





This work 


Kiselev [4] a 


IKP [5] 


WSL [27] 


DSV [28] 


DW [29] h 


pB c ^D 


o u +0 - 01 

U - 1 ^-0.02 


0.32 [0.29] 


0.189 


0.16 


0.075 


0.255 


pB c ^D s 


19+0- 02 

u - iy -0.01 


0.45 [0.43] 


0.194 


0.28 


0.15 


- 


A B C -+D* 


12+ - 02 

u - iz -0.01 


0.35 [0.37] 


0.133 


0.09 


0.081 


0.257 


A B ^ D *s 
A 


17+0- 01 

u - i '-0.01 


0.47 [0.52] 


0.142 


0.17 


0.16 


- 



a The non-bracket (bracketed) results are evaluated in sum rules (potential model) 
6 We quote the result with ui = 0.7GeV 



Table HIl The error is from the combined uncertainty in the hadronic parameters: (1) the 
shape parameters: u>b = 0.60 ±0.05 for B c meson wave function, dp = (0.5±0.1)GeV for 
D(*> meson and ao a = (0.4 ± 0.1)GeV for D s * meson wave function |20|; (2) the decay 
constants in the wave functions of charmed mesons, which are given in Table HI Since the 
uncertainties from decay constants of D( s ) and the shape parameters of the wave functions 
are very small, the relevant uncertainties to the form factors are also very small. We can 
see that the SU(3) symmetry breaking effects between B c to Z)W and B c to D s * form 
factors are large, as the decay constant of D s is about one-fifth larger than that of the D 
meson. 

n the literature there are already lots of studies on B c — > D\*l transition form factors 



(•) 
5|, |27f|29|, whose results are collected in Table ILT1 Our results are generally close to 

the covariant light-front quark model results of [27] and the constituent quark model 

results of [5|. However, other results collected in Table HH especially for the QCD sum 

rules (QCDSR) [4| and the Bauer, Stech and Wirbel (BSW) model I2J deviate a lot 



10 



numerically. The predictions of QCDSR [4J are larger than those in other works 



[4J are larger than those in other works [5|, 1271 



29|. The reason is that they have taken into account the a s /v corrections and the form 



factors are enhanced by 3 times due to the Coulomb renormalization of the quark-meson 



vertex for the heavy quarkonium B c . The results of BSW model [28j are quite small due 
to the less overlap of the initial and final states wave functions. Although, the included 
flavor dependence of the average transverse quark momentum in the mesons can enhance 
the form factors for B c — > D? s <, transitions, their predictions are still smaller than other 
models. The large differences in different models can be discriminated by the future LHC 
experiments. 

B. Branching Ratios 

With the decays amplitudes A obtained in SecJTTl the branching ratio BIZ reads as 

Bn = §St B ^ {l - (ra + ra)2)(1 - (r2 - r3)2)| - 4 ' 2 - (16) 

As stated in Sec [Til the contributions from the penguin operators are absent, since the 
penguins add an even number of charmed quarks, while there is already one from the initial 
state. There should be no CP violation in these processes. We tabulate the branching 
ratios of the considered decays in Table UTTI and IIVI The processes (l)-(4) in Table UTTI have 
a comparatively large branching ratios (10~ 5 ) with the CKM factor V* b V u d ~ A 2 . While the 
branching ratios of other processes are relatively small due to the CKM factor suppression. 
Especially for the processes (l)-(4) in Table HVl these channels are suppressed by CKM 
element V u b/V c b and V c d/V u d- Thus their branching ratios are three order magnitudes 
smaller. 

For comparison, we also cite other theoretical results |4j, |5|, |7|, |8|, [l0| for the double 



charm decays of B c meson in Tables IIHI and IIVI In general, the results of the various 
model calculations are of the same order of magnitude for most channels. However the 
difference between different model calculations is quite large. This is expected from the 
large difference of input parameters, especially the large difference of form factors shown 
in Table HTJ As stated in the introduction, all the calculations of these B c to two D 
meson decays in the literature use the same naive factorization approach. Their difference 
relies only on the input form factors and decay constants. Therefore the comparison of 
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TABLE III: Branching ratios (10~ 6 ) of the CKM favored decays with both emission and an- 
nihilation contributions, together with results from other models. The errors for these entries 
correspond to the uncertainties in the input hadronic quantities, from the CKM matrix elements, 
and the scale dependence, respectively. 





channels 


This work 


Kiselev[4 


IKP[5' 


IKS [7] 


LC[8l 


CFflO] 


1 


B c -> D+D° 


on+6+1+2 
OZ -6-l-4 


53 


32 


33 


86 


8.4 


2 


B c ->• D+D*° 


0,, +7+2+3 
04l -6-l-3 


75 


83 


38 


75 


7.5 


3 


B c -»• D*+D° 


-1 9 +3+l+0 
iz -3-0-l 


49 


17 


9 


30 


84 


4 


B c -> D*+D*° 


o. +9+2+0 
-J 4 -8-l-0 


330 


84 


21 


55 


140 


5 


B c -» D+L>° 


9 q+0.4+0. 1+0.2 
z -°-0.4-0.1-0.2 


4.8 


1.7 


2.1 


4.6 


0.6 


6 


5 C ->• D+D*° 


9 c+0.4+0. 1+0.1 
z -°-0.6-0. 1-0.2 


7.1 


4.3 


2.4 


3.9 


0.53 


7 


5 C ->• D*+5° 


n 7+O.I+O.O+O.O 
u -' -0.2-0.0-0.0 


4.5 


0.95 


0.65 


1.8 


5 


8 


5 C -> D*+D*° 


9 O+0.7+0. 1+0.1 
z -°-0.6-0. 1-0.0 


26 


4.7 


1.6 


3.5 


8.4 



TABLE IV: Branching ratios (10 -7 ) of the CKM suppressed decays with pure emission contri- 
butions, together with results from other models. The errors for these entries correspond to the 
uncertainties in the input hadronic quantities, from the CKM matrix elements, and the scale 
dependence, respectively. 





channels 


This work 


Kiselev[4 


IKP[5 


IKS [7] 


1 


B c 


^D+D° 


-, n+o.2+0. l+o.o 
■•••"-o.i-o.o-o.o 


3.2 


1.1 


3.1 


2 


B c 


-> D + D*° 


n 7+0.1+0.1+0.0 
u -' -0.2-0.0-0.0 


2.8 


0.25 


0.52 


3 


B c 


-+ D*+D° 


n Q+0.1+0. 1+0.0 
u - J -0.2-0. 0-0.0 


4.0 


3.8 


4.4 


4 


B c 


v r)*+ r)*0 


n n+0.2+0. 1+0.2 
u -°-0.1-0. 0-0.0 


15.9 


2.8 


2.0 


5 


B c 


-► D+D° 


o n +5+3+l 
ou -4-2-l 


66 


25 


74 


6 


B c 


-»• D+D*° 


1 q+3+2+0 
iy -3-l-l 


63 


6 


13 


7 


B c 


-^ D* s +D° 


9C-+4+2+0 
zd -3-2-l 


85 


69 


93 


8 


B c 


-»• D*+D*° 


9/1+3+2+1 
z4 -3-2-l 


404 


54 


45 



results with any of them is straightforward. Larger branching ratios come always with 
the larger form factors. As stated in the previous subsection, our results of form factors 
are comparable with the relativistic constituent quark model (RCQM) p, [D], thus our 
branching ratios in Table ILTT1 are also comparable with theirs except for the processes B c — > 
D*+D*° and B c — > D* + D*°. Due to the sizable contributions of transverse polarization 
amplitudes, our branching ratios are larger than those in RCQM model, whose transverse 
contribution is negligible. 
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Since all the previous calculations in the literature are model calculations, it is dif- 
ficult for them to give the theoretical error estimations. In our pQCD approach, the 
factorization holds at the leading order expansion of m£>/mB- At this order, we can 
do the systematical calculation, so as to the error estimations in the tables. The first 
error in these entries is estimated from the hadronic parameters: (1) the shape pa- 
rameters: lob = 0.60 ± 0.05 for B c meson, ao = (0.5 ± 0.1)GeV for D^-*' meson and 



Q-Ds = (0.4±0.1)GeV for D s * meson 20j; (2) the decay constants in the wave functions of 
charmed mesons, which are given in Table HI The second error is from the uncertainty in 
the CKM matrix elements, which are also given in Table [I] The third error arises from the 
hard scale t varying from 0.75t to 1.25t, which characterizing the size of next-to-leading 
order QCD contributions. The not large errors of this type indicate that our perturba- 
tive expansion indeed hold. It is easy to see that the most important uncertainty in our 
approach comes from the hadronic parameters. The total theoretical error is in general 
around 10% to 30% in size. 

The eight CKM favored channels (proportional to \V c b\) in Table HTT1 receive contribu- 
tions from both emission diagrams and annihilation diagrams. From FigJTJ one can find 
that the contributions from the factorizable emission diagrams are color-suppressed. The 
naive factorization approach can not give reliable predictions due to large non-factorizable 



contributions 30] . As was pointed out in SecJTTl the non-factorizable emission diagrams 
give large contributions in pQCD approach because the asymmetry of the two quarks 
in charmed mesons. Thus, the branching ratios of these decays are dominated by the 
non-factorizable emission diagrams. 

The eight CKM suppressed channels (proportional to \V u b\) in Table HVl can occur only 
via emission type diagrams. There are two types of emission diagrams in these decays, one 
is color-suppressed, one is color favored. It is expected that the color-favored factorizable 
amplitude T e z dominates in eq.f TTS]) . However, the non-factorizable contribution A4 e 2, 
proportional to the large C*2, is enhanced by the Wilson coefficient. Numerically it is 
indeed comparable to the color- favored factorizable amplitude. This large non-factorizable 
contribution has already been shown in the similar B — > Dtt decays theoretically and 
experimentally [25]. In all of these channels the non-factorizable contributions play a 
very important role, therefore the branching ratios predicted in table 



not like the previous naive factorization approach calculations 
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a table I 

[ y, y, k 



H and lIVl are 
10| . They are 



TABLE V: The transverse polarizations fractions (%) for B c — > VV. The errors correspond to 
the uncertainties in the hadronic parameters and the scale dependence, respectively. 



K T 



B r -> D* + D*° 



58 -3-0 



B r ->• D*+D 



*0 



68+n 



B r -> D*+D*° 



4±i±i 



s r -> £>:+l> 



*o 



6 +l+2 
°-0-l 



not simply proportional to the corresponding form factors any more, but with a very 
complicated manner, since we have also additional annihilation type contributions. 

From Table III and IV, one can see that as it was expected the magnitudes of the 
branching ratios of the decays B c — , D+D° and B c — > DfD° are very close to each 
other. In our numerical results, the ratio of the two decay widths is estimated as 
J- c ^ n+noi ~ 1-3- They are very suitable for extracting the CKM angle 7 though the 
amplitude relations. Hopefully they will be measured in the experiments soon. How- 
ever, the decays B c — > D + D°, D + D° are problematic from the methodic point of view for 
BK(B C -* D + D°) < BTZ(B C -* D + D°). The corresponding ratio in B c -* D + D°,D + D° 
decays is r \ B c ^ D +i)0) ~ 1CT 3 , which confirm the latter decay modes are not useful to 
determine the angle 7 experimentally. 

For the B c decays to two vector mesons, the decays amplitudes A are defined in the 
helicity basis 

A= J2 iai 2 , ( 17 ) 

i=0,+ - 

where the helicity amplitudes Ai have the following relationships with A L ' N,T 

A = A L , A ± = A N ±A T . (18) 

We also calculate the transverse polarization fractions IZt of the B c — > DJsD* decays, 
with the definition given by 

\AJ 2 + \A I 2 

n T = — — — fi9) 

T \Ao\* + \A+\ 2 + \A-\ 2 ' [ } 

This should be the first time theoretical predictions in the literature, which are absent 
in all the naive factorization calculations. According to the power counting rules in 
the factorization assumption, the longitudinal polarization should be dominant due to 
the quark helicity analysis. Our predictions for the transverse polarization fractions of 
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the decays B c — > DJtD* , which are given in Table |V] are indeed small, since the two 
transverse amplitudes are down by a power of t2 or r^ comparing with the longitudinal 
amplitudes. However, for B c —¥ D*AD*° decays, the most important contributions for 
these two decay channels are from the non-factorizable tree diagrams in Fig. [TJc) and[Tfd). 
With an additional gluon, the transverse polarization in the non-factorizable diagrams 
does not encounter helicity flip suppression. The transverse polarization is at the same 
order as longitudinal polarization. Therefore, we can expect the transverse polarizations 
take a larger ratio in the branching ratios, which can reach ~ 60%. The fact that the 
non-factorizable contribution can give large transverse polarization contribution is also 
observed in the B° — > p°p°, ww decays [3l| and in the B c — > D* + u decay |32| . 

IV. CONCLUSION 

All the previous calculations in the literature for the B c meson decays to two charmed 
mesons are based on the very simple naive factorization approach. The branching ratios 
predicted in this kind of model calculation depend heavily on the input form factors. 
Since all of these modes contain dominant or large contributions from color-suppressed 
diagrams, the predicted branching ratios are also not stable due to the large unknown 
non-factorizable contributions. In this paper, we have performed a systematic analysis of 
the double charm decays of the B c meson in the pQCD approach based on kx factoriza- 
tion theorem, which is free of end-point singularities. All topologies of decay amplitudes 
are calculable in the same framework, including the non-factorizable one and annihilation 
type. It is found that the non-factorizable emission diagrams give a remarkable contri- 
bution. There is no CP violation for all these decays within the standard model, since 
there are only tree operators contributions. The predicted branching ratios range from 
very small numbers of (9(10~ 8 ) up to the largest branching fraction of (9(10 -5 ). Since all 
of the previous naive factorization calculations did not give the theoretical uncertainty 
in the numerical results, it is not easy to compare our results with theirs. The theo- 
retical uncertainty study in the pQCD approach shows that our numerical results are 
reliable, which may be tested in the upcoming experimental measurements. We predict 
the transverse polarization fractions of the B c decays with two vector D* mesons in the 
final states for the first time. Due to the cancelation of some hadronic parameters in 
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the ratio, the polarization fractions are predicted with less theoretical uncertainty. The 
transverse polarization fractions are large in some channels, which mainly come from the 
non-factorizable emission diagrams. 
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Appendix A: Factorization formulas for B c — > VV 

In the B c decays to two vector meson final states, we use the superscript L, N and 
T to denote the contributions from longitudinal polarization, normal polarization and 
transverse polarization, respectively. For the CKM favored B c — > DJ^D* decays, the 
decay amplitudes for different polarizations are 

/2 f 1 f°° b 2 cu 2 

^C F f B f 3 nM 4 B l dx 2 l b 1 b 2 db 1 db 2 M^)exp(-^)x 

{[-(r 2 - 2)r b + 2r 2 x 2 - x 2 ]a s (t a )h e (a e , (3 a , h, b 2 )S t (x 2 ) exp[-S ab (t a )] 

+rla s (t b )h e (a e , (3 b , b 2 , 6i)S' t (a:i) exp[-S ab (t b )]}, (Al) 



b u> 

b 2 b 3 db 2 db 3 <p 2 (x 2 )(p 3 (x 3 ) exp( — 2 — -) x 



M L e = --CpfBTiM 4 / dx 2 dx 3 \ 
<J Jo Jo 

{[1 - xi - x 3 + r 2 (l - x 2 )]a s (t c )h e ((3 c , a e , b 3 , b 2 ) exp[-S cd (t c )] - 

[1 - xi - x 2 + x 3 - r 2 (l - x 2 )]a s (td)h e (Pd, ®e, h, h) exp[-S' cd (^)]}, (A2) 

/»1 POO 

J^a = -^C F fB^M B I dx 2 dx 3 I b 2 b 3 db 2 db 3 (j) 2 (x 2 )(j) 3 (x 3 ) x 

Jo Jo 

{[1 - x 2 }a s (t e )h e (a a , f3 e , b 2 , b 3 ) exp[-S ef (t e )]S t (x 3 ) - 

[1 - x 3 ]a s {t f )h e (a a , (3 f , b 3 , b 2 ) exp[-S ef (tf)]S t (x 2 )}, (A3) 
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h 2 , ! 2 



'B< 



-pN __ _ 



o pi poo 

a = oCf/b-kMs / dx 2 dx 3 / b 1 b 2 db 1 db 2 (/)2(x2)(f>3(x3) exp 
•J Jo Jo 

{[xi +x 3 -l - r c ]a s (t g )h e (f3g,a a ,b 1 ,b 2 ) exp[-S gh (t g )] 

-[r b - x 2 ]a 8 (th)h e (p h ,aa,b 1 ,b 2 ) exp[-S gh (t h )]}, (A4) 

^C F f B f 3 r 3 irM B J dx 2 J b 1 b 2 db 1 db 2( j )2 (x 2 )exp(-^)x 
{[2 - r b + r 2 (4r 6 - x 2 - l)]a s (i a )/i e (a e , /3 a , 61, b 2 )S t (x 2 ) exp[-S ab (t a )] 
-r 2 a s (t b )h e (a e , (3 b , b 2 , b^S^x^ exp[-S ab (t b )]}, (A5) 

^C F f B f 3 r 3 nM B I dx 2 l b 1 b 2 db 1 db 2( f> 2 (x 2 )exp(-^-)x 
{[2 - r b - r 2 (l - x 2 )]a s (£ a )/i e (a e , &, 61, &2)-S't(a; 2 ) exp [-S^a)] 
-r 2 a s (t b )h e (a e , (3 b , b 2 , bjStfa) exp[-5' a6 (^)]}, (A6) 

8 C 1 f°° h 2 u 2 

Mf = -M T e =-C F f B nM A B r 3 / dx 2 dx 3 / b 2 b 3 db 2 db 3 <f> 2 (x2)cf>3(x3)eM-^7r-) 
3 Jo Jo * 

x{[xi +x 3 - l]a s (t c )/i e (/3 c ,a e ,6 3 ,6 2 )exp[-S' crf (t c )] - 

[xi - x 3 ]a s (t d )/i e (/3 d , a e , 63, 62) exp[-S' cd (i d )]}, (A7) 

-8C F f B nM B r 2 r 3 / dx 2 dx 3 / b 2 b 3 db 2 db 3 (/) 2 (x 2 )(j) 3 (x 3 ) x 

Jo Jo 

{[2 - :r 2 ]a s (t e )/i e (a a , /3 e , 6 2 , 63) expf-^/^)]^^) - 

[2 - x 3 ]a s (i/)/i e (a a , /?/, 63, 6 2 ) exp[-S , e/ (t / )]S , t (a; 2 )}, (A8) 

J"J = -8C F f B irM B r 2 r 3 / <ir 2 <ir3 / b 2 b 3 db 2 db 3 (f)2(x 2 )(f) 3 (x 3 ) x 

Jo Jo 

{x 2 a s (t e )/i e (a a , /3 e , 62, 63) exp[-S ef (t e )]S t (x 3 ) + 

x 3 a s (t f )h e (a a , f3 f , b 3 , b 2 ) exp[-S ef (tf)]S t (x 2 )}, (A9) 

8 r 1 C°° b 2 u] 2 

■Ma = nC F f B 7rM B / dx 2 dx 3 / bifodbidhifaMfofa) exp( — 1 —^-) x 
<J Jo Jo ^ 

{[r 2 2 {x 2 - 1) + r 2 (rr 3 - l)]a s (i g )/i e (^ g , a a , b ± , b 2 ) exp[-S gh (t g )] 

-\r\x 2 + r\x 3 - 2r 2 r 3 r b ]a s (t h )h e ((3 h , a a , 61, 6 2 ) exp[-S'^(^)]}, (A10) 

-Ml = 7>C F f B nM B I dx 2 dx 3 / 6 1 6 2 rf6 1 rf6 2 2 (a; 2 )03(a;3)exp( — y^-) x 
{[r 2 (> 2 - 1) -rj(x 3 - l)]a s (t g )h e (/3 g , a a ,b!,b 2 ) exp[-S gh (t g )] 
-\r\x 2 - rlx 3 ]a s (t h )h e (f3 h ,a a ,biM) ex-V[-Sgh{th)]}- (All) 



-pN 
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For the CKM suppressed B c — > DJtD* decays, the decay amplitudes for different 
polarizations are 



-C F f B f 3 7tM B j dx 2 J b 1 b 2 db 1 db 2 c} )2 (x 2 )eM- A ^-: 



{[-(r 2 - 2)r b + 2r 2 x 2 - x 2 ]a s (t a )h e (a e , (3 a , b u b 2 )S t (x 2 ) exp[-S ab (t a )] 
+r 2 2 a s {t h )h e {a e , (3 b , b 2 , b^S^x^ exp[-Sab(tb)]}, (A12) 



T. 



/2 f 1 f°° b 2 uj 2 

^C F f B f 3 r 3 7rM B / dx 2 / 6 1 6 2 rf6 1 d6 2< /) 2 (a; 2 ) exp( — 1 —^-) x 

{[2 - r b + r 2 (4r 6 - x 2 - l)]a s (t a )h e (a e , (3 a , b t , b 2 )S t {x 2 ) exp[-S ab (t a )] 

-r 2 a s (t b )h e (a e , (3 b , b 2 , b^S^x^ exp[-S ab (t b )]}, (A13) 

/2 Z" 1 /"°° b 2 cu 2 

^C F f B f 3 r 3 7rM B J dx 2 j 6i6 2 d6 1 d6 2 ^ 2 (x 2 ) exp(-^) x 

{[2 - r 6 - r 2 (l - x 2 )]a s (t a )/i e (a e , /3 a , 61, b 2 )S t (x 2 ) exp[-S ab (t a )] 

~r 2 a s (t b )h e (a e , f3 b ,b 2 ,b 1 )S t (x 1 ) exp[-S ab (t b )]}, (A14) 



Q /"l /"OO 

A4^ 2 = -C F f B rcM B / rfx 2 (ia;3 / b 2 b 3 db 2 db 3 (f) 2 (x 2 )(f) 3 (x 3 ) exp( 
<J Jo Jo 



7l 2 / ! 2 



X 



2 
{[2 - zi - £2 - x 3 - r 2(! - z 2 )]a s (t c )/i e (/3 c , a e , 6 3 , b 2 ) exp[-S cd (t c )] - 

[x 3 - xi + r 2 (l - z 2 )]a s (t d )/i e (/3 d , a e , 6 3 , 6 2 ) exp[-5' cd (^)]}, (A15) 



°1 W B' 



Q pi pOO 

M%, = ~Ml 2 = -C F f B ixM B \ dx 2 dx 3 I b 2 b 3 db 2 db 3 (j) 2 (x 2 )(f) 3 (x 3 ) exp( 
-J Jo Jo 

x{M^i -z 3 )]a s (t c )/i e (/3 c ,a e ,6 3 ,6 2 )exp[-S' cd (t c )] + 

[2r c - r 3 (l - xi - Xs)]a«(*d)MAi» a e , b 3 , b 2 ) exp[-S cd (t d )]}. (Af 6) 

Appendix B: Scales and related functions in hard kernel 

We show here the functions h e , coming from the Fourier transform of hard kernel, 

h e (a,PMM) = h 1 (a,bi)xh 2 (P,b 1 ,b 2 ), 

h x (aM) = /^(v^i). a>0 

[ K (iy/—abi), a < 

h (Rb b\ = i 6{hl ~ MMvW^vW + (61 <* 62), /? > m) 



where Jo is the Bessel function and Kq, Iq are modified Bessel function with Ko(ix) = 
^(—Nq(x) +iJo(x)). The hard scale t is chosen as the maximum virtuality of the internal 
momentum transition in the hard amplitudes, including l/bi(i = 1,2,3): 

t a = max(Vl(^|,\/pJ,l/&i,l/&2), h = max(A/ja~[, ^/\^\,l/b 1 ,l/b 2 ), 

t c = max(v / faT[, \/pc[> V & 2, 1/63), U = max(v / M, vJA*I> I/62, I/63), 

t e = max(\/|c^[, VW, I/62, l/M) ^/ = m ax(A/|aJ, y 1/3/1,1/62,1/63), 

t 5 = max(\/[aj, y |)9 5 |,l/6i,l/62), t fc = max(A/jo~|, a/I/^L I/61, l/6 2 ), (B2) 

where 

a e = (1 - a: 2 )(xi - r 2 2 )(l - r 3 2 )il4, a a = -(1 + (r 3 2 - l)x 2 )(l + (r 2 2 - l)a; 3 )M 2 , 

Pa = VI + (r 2 - l)(x 2 + r 2 (l - x 2 ))]M 2 , ft = (1 - r 2 )^ - ri)A*l, 

/?c = [r c 2 - (1 - x 2 (l - r 2 ))(l - Xl - x 3 (l - r 2 ))]]M 2 , 

/?d = (1 - x 2 )(l - r 2 )[x! - x 3 - r 2 (l - x 3 )}M 2 B , 

ft = -[1 + (r 2 - l)x 2 ]M 2 , 0, = -[1 + (r 2 - l)x 3 ]M 2 , 

ft = [r 2 + (1 - x 2 (l - r 2 ))(x! + x 3 - 1 - r 2 x 3 )]M 2 , 

& = [r 2 - x 2 (r 2 - l)( Xl - x 3 (l - rl))]Ml. (B3) 

The Sudakov factors used in the text are defined by 



S ab {t) --■ s(^£xi,6i) + s(^£x 2 ,6 2 ) + - / — 7,(/x) + 2 / — 7«(a0> 

6 Jl/bi V Jl/b 2 V 



.Mb , , , ,M B 



M B M B M B 

S'cd(t) = s(— =0:1,62) + s(— ; =x 2 M) + s(—, =x 3 ,b 3 ) 

6 Jl/bi f 1 Jl/bs ^ 

Sef(t) = s(^X 2 ,b 2 ) + s(^xM+2 f ^7^+2/* % q (jM), 
V2 V2 Ji/b 2 ^ Jl/h ^ 



S 9 h{t) = s(—^x 1 ,b 1 ) + s{—^X2,b2)+s(—^x 3 ,b 2 ), 



■ M B M b M b 

-J^Xl, 61) + S \—fr X ^ & 2) + S{— = 

5 f l dfi /"* d/j, 



„ 7 9 (/i)+4 / — 7(? (/i), (B4) 

■J 7 1/61 A 1 J 1/62 f 1 



where the functions s(Q, 6) are defined in Appendix A of 24J. 7^ = —a s /iT is the anoma- 
lous dimension of the quark. 
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Appendix C: Meson Wave functions 



In the nonrelativistic limit, the B c meson wave function can be written as 



®rJx) 



JJb 



[GP + M Bc ) 75 (5(x-r c )]exp(- 



2, ,2 



b 2 LO 



(CI) 



in which the last exponent term represents the k^ distribution. Here, we only consider 
the dominant Lorentz structure and neglect another contribution in our calculation 34]. 
In the heavy quark limit, the two-particle light-cone distribution amplitudes of 
D( s )/D?s meson are defined as [35 1 

• i 



(D {s) (P 2 )\q a (z)c^0)\0) 
(D{ s) (P 2 )\q a (z)c (O)\O) 



V2N C 



c JO 



V2N< 



c JO 



dxe lxy2 - z fr 5 (P 2 + m DM )<p D( Jx,b)]ap, 



V s ) (s) 



+fr( p 2 + m D * )<f) T D *(x, b)] a p, 

\ s ) (s) 



(C2) 



with the normalization conditions: 



dx<& u {x,0)= / n 1 dx^ f . ) (x,0) = ^ 



(C3) 



where we have assumed fo* = fn* • Note that equations of motion do not relate </>£ 



W 



(«) 



and (fi)* ■ We use the following relations derived from HQET J36( to determine /d 



w 



oo 



(«) 



./; 



'% 



D 



(«) 



M 



m D 



fo {s) ■ 



(C4) 



M 



The distribution amplitude </> A, is taken as [l 



i? 



(*) 
w 






f D wx(l -x)[l + a D w(l -2z)]exp(- 



h 2 i i 2 



(») 



)• 



(C5) 



'« 2 

We use a D = 0.5 ± 0.1, co D = O.lGeV for D/D* meson and a D = 0.4 ± 0.1, u Da = 0.2GeV 
for D s /D* meson, which are determined in Ref. 20] by fitting. 
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